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Abstract 

Exclusive photoproduction of heavy quarkonia at HERA is analyzed in the frame- 
work of the Regge-eikonal approach together with the nonrelativistic bound state 
formalism. Total and differential cross-sections for the process 7 +p — > (QQ)is +P 
are calculated. The model predicts cross-sections of Exclusive Double Diffractive 
Events (EDDE) at TeVatron and LHC. 
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1 Introduction 



The study of properties of bound states of heavy quarks plays a central role in the under- 
standing of strong interactions and verification of different QCD inspired and nonpertur- 
bative models, because such processes give a very exciting possibility to observe interplay 
of "hard" and "soft" regimes [1]. There have been intensive experimental studies of the 
J/ty and T production in ep collisions at HERA [2], [3] and also a lot of theoretical inves- 
tigations [4]- [8]. 

In this paper we consider exclusive photoproduction of V = QQis states from another 
viewpoint. There are some other interesting processes that will be investigated at present 
and future hadronic colliders. We need estimations of cross-sections for such processes and 
can use the data from HERA as a source of normalization of phenomenological models. 
Here we show that the extended Regge-Eikonal approach [9]- [12] gives not only a good 
description of the data on exclusive vector meson photoproduction but can be used also 
to predict rates of Exclusive Double Diffractive Events (EDDE) at LHC and Tevatron. 
The advantages of these events have been considered in [13], [14]. 



2 Calculations 

In Fig. 1 we illustrate in detail the process 7(g) +p(p) — > V{p v ) +p(p'). Off-shell proton- 
gluon amplitude T in Fig. 1 is treated by the method developed in Ref. [9], which is based 
on the extension of Regge-eikonal approach, and succesfully used for the description of 
the data from hadron colliders [10]-[12]. The amplitude A of the process 7(g) + g(Ki) — > 
V(p v ) + g(n 2 ) is calculated in the nonrelativistic bound state approximation(see [4]-[6] 
and ref. therein): 
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where pi = M%, cq is the charge of heavy quark, R v0 is the absolute value of the vector 
meson radial wave function at the origin, e U)7 are photon and vector meson polarization 
vectors correspondingly. Permutations are taken for all gauge bosons. Notations and 
vector decompositions that are used in the article are 



«i = K + — , k 2 = k- — , p = p + -^q , q = q - —p , (3) 



Q = -q , p ' —rrip , q =p =0 , s~2pg , 
K = Tr( a P + /3q) + n± , x v = , k ± = -k , y = —r , y = 



2 v y ' f HJ ' " ' v s ' x ' y Ml ' y " Ml ' 

O 2 - A 2 

A = x v ([1 + y Q + y A ]p' - y A q') + A ± , y Q = , t ~ A 2 ± = -A 2 , y A = — , 
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p v = q + A,y = — 

V 

Photon and vector meson polarization vectors in the general case (Q ^ 0) can be repre- 
sented as follows: 

e 7± g = e 70 g = , e 7 ^_ = -e 7 2 = -1 , e 70 = — (q' + x v y Q p') , (4) 

2(vA) . . 2 _ 2 

e v ±Pv = e v \\Pv = , e v± = v ± H (p - q ) , v ± = -u , 

s 

= jfW - x v (l - Va)p + Aj_) 
For the amplitude of the process 7(<2 r ) — > V(p„) we have: 

M -/(^ ( K ; + ,:o)( K i + i o) - 4 °"°^- < 5 » 

— 9af3 , (/ ) 

KiK 2 

a = p «2 , "2 = v K i , (8) 

Generally the amplitude T® p ^ gp can be represented in the Regge-eikonal form [10], [12] 
with fixed parameters of trajectories from Ref. [12] (see Table. 1), in which the eikonal is 
dominated by three vacuum trajectories (Pomerons with different properties). It follows 
from the analysis below that at small t the amplitude T®^ gp takes the simple Regge form, 
which is dominated by the 3rd ("hard") Pomeron: 

T n „ c ( e ^J^_X P ' A(t) e 6 ° 3)< bf = ^ + °' 5r ^ (9) 

1 gp->gp — C 9P Sq _ k 2J e ' W ^ > 

where s ^ 1 GeV is the scale parameter of the model that is used in the global fitting 
of the data on pp{pp) scattering [11], [12], r pPs , ap 3 (t) = op 3 (0) + a'p 3 (0)t are defined in 
Table. 1, r 2 gPi and c gp are extracted by the procedure (16)-(21). With notations (3) we 
have: 

M 4 x v M 4 x v da 

d k — 7r — ^—^-dadfldy = —71 — - — - — dy'dy (10) 
32 64 a 

In the limit Q — > , t — > only the amplitude M±± survives: 



/ c \ 2a P (0) 

M ± J 2 * K2 Mt) 2 4 P {£) e 2 ^, (11) 

h = &!, 3) +c4(0)m- (12) 

So 
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Now let us extract the values of parameters from the fit to the data on elastic J/^f 
photoproduction [2]. At first we write the amplitude M±± in the Regge-eikonal form 
with parameters from Table. 1 and the coefficient that corresponds to the simple Vector 
Dominance Model (VDM): 



M ±± = \ { J db 2 J (bV^tf- -, (16) 

V OL e M v Jo 1% 



where 



JS) / e \ap.(0) — 1 p pf 

S^&U-**!) e — 2 , (17) 

s V s J 4tt pf 



p 2 = 4c^(0) In (e"**^) + r 2 gPi + 0.5r 2 Pi 

As will be seen below, in our case the VDM plus Regge-eikonal approach representa- 
tion (16) is applicable. 

Results of this fit for J/^> meson are shown in Figs. 2-5. As we see from figures, the 
main contribution to the cross-section is given by the Born term of the 3rd Pomeron. The 
1st "soft" Pomeron gives no contribution. The term corresponding to the 2nd Pomeron 
vanishes faster with t, and gives the contribution less than 1%, when t < —0.2 GeV 2 . 
Numerical estimations show that absorbtive corrections play minor role at t ~ t* — 
= — 1/26 3 , where 63 is obtained from (12). Using these facts, we keep in (16) only the 
Born term for the 3rd Pomeron with parameters 

r 2 gP3 = 2.54 ± 0.41GeV" 2 , Cj % p = 1.11 ± 0.07 , X 2 /dof = 1.48 (18) 

and take the integral at t — t*. Now we can estimate the constant c gp in (11) from the 
comparison of two formulae for the amplitude M^±: 



V a e M v (3) d7r r (3) ( 1 n\ 

K v I v (t*) °<* Z2a s I v (t*)VK^3 vp ' 1 } 

Taking for J/ty mesons 

Mji^ = 3.1 GeV , a s (M 2 ^) = 0.25 , (20) 

Ij/v(t*) ~ 0.83 , 35 GeV < W = V* < 260 GeV , 

r(J/tf -> e + e~) = 5.26 ± 0.37 keV , K NLO ~ 2 (see, for example, [15]) 



4 



we get from (19): 



c gp = 3.5 ± 0.4 



(21) 



Here errors are estimated from uncertanties of quantities in (19). 

The data on T production [3] gives the possibility to check the model predictions. The 
result of ZEUS collaboration for the ratio of total cross-sections of J/ty and T photopro- 
duction: 



a ^ Tp - (4.8 ± 2.2(stat.)+^(sys.)) • 10~ 3 



(22) 



Gyp— ► 

If we assume that the constant c gp is the same for both processes, and the slope of the 
exponent does not change much with energy, then from the expression (11) we will get at 
the same value of W: 



Gyp^Tp 
Gyp-^J/Wp 
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(23) 



T(T -> e + e~) = 1.32 ± 0.04 ± 0.03 keV , 
My = 9.46 GeV , a s (M r ) ~ 0.2 , I r ~ 0.21 , 

k nlo ~ z — (see Ref. [15]) , 



(24) 
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and uncertainty of the result originates from the errors of parameters in (23). Theoretical 
estimation does not contradict the experimental value (22). 

The second estimation can be done for the EDD dijet production at TeVatron energies. 
Recent CDF results [16], [17] for the upper bound of the cross-section of the process p+p — > 
p + jet + jet + p are the following: 



E T >7 GeV , a < 3.7 nb , (25) 
E T > 10 GeV , a < 0.97 ± 0.065 (stat.) ± 0.272(sys.) nb , 
E T > 25 GeV , a < 34 ± 5 (stat.) ± 10(sys.) pb . 

After theoretical calculations by the method developed in Refs. [13], [18] we extract upper 
bounds for the parameter c gp from (25): 

E T >7 GeV , c gp < 3.3 , (26) 
E T > 10 GeV , c gp < 3.4 , 
E T > 25 GeV , c gp < 4.2 . 

Values of c gp are close to our estimation (21). 
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Conclusions 

We can conclude that the generalized Regge-eikonal approach with 3 different Pomerons 
describes well the data on J/^> production. The main contribution to the cross-section 
comes from the term corresponding to the 3rd, so called, "hard" Pomeron. This makes 
possible to extract the corresponding parameter of the model. 

The upper bound for the same parameter is found to be close to our result, when 
calculated from experimental estimations on EDD di-jet production made by CDF. It 
indicates once more the applicability of the Regge-eikonal approach and gives us the tool 
for further predictions. 
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Table 1.: Parameters 0^(0), a' P .(0), r 2 Pi are obtained from the fit to the data on 
pip) + P ~ * Pip) +P [12] and remain fixed during the J/\l/ data fitting. 



Pomeronj 


1 


2 


3 


a Pi (0)-l 


0.0578 ± 0.0020 


0.1669 ±0.0012 


0.2032 ±0.0041 


o'p.(O) (GeV- 2 ) 


0.5596 ± 0.0078 


0.2733 ±0.0056 


0.0937 ±0.0029 


r 2 pPi (GeV- 2 ) 


6.3096 ± 0.2522 


3. 1097± 0.1817 


2.4771 ± 0.0964 
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Figure captions 

Fig. 1: Diagram for the process 7 + p — > V + p. 
Fig. 2: Diagrams for the process 7 + g* — > V + g*. 

Fig. 3-5: Differential cross-sections of the process 'j+p^V+p&t different values of W. 
Solid line is the Born term for the 3rd Pomeron and dashed curve is the unitarized 
result. 



Fig. 6: Total cross-section of the process 7 + p — > V + p. Solid line is the Born term for 
the 3rd Pomeron and dashed curve is the unitarized result. 
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Figure 4: 
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Figure 5: 
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